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_C ■ Abstract. The spectrum of the Hermitian Hamiltonian H = p 2 + V(x) is real and 

discrete if the potential V(x) — > oo as x — > ±oo. However, if V(x) is complex and VT- 
symmetric, it is conjectured that, except in rare special cases, V(x) must be analytic in 
order to have a real spectrum. This conjecture is demonstrated by using the potential 

pi \ V(x) — (ix) a \x\ h , where a, b are real. 
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q ■ 1. Introduction 

The field of VT quantum mechanics [1] has attracted significant interest in recent years 
and a large community of active researchers has developed. This area of study began 
with the observation that the complex PT-symmetric non-Hermitian Hamiltonian 



ctf 



><■ H = v 2 + x 2 (ix) e (e>0) (1) 

has a positive real discrete eigenspectrum [2j [3] . The property of VT symmetry is not 
sufficient to guarantee that the eigenvalues of a non-Hermitian Hamiltonian are real; VT 
symmetry merely ensures that the secular determinant det(H — IE) is a real function 
of E [5]. The eigenvalues of H are the roots of 

det(# - IE) = 0, (2) 

and thus the condition of VT symmetry implies that the eigenvalues are either real or 
come in complex-conjugate pairs. If the eigenvalues of a PT-symmetric Hamiltonian 
are all real, we say that the Hamiltonian has an unbroken VT symmetry, but if there 
are any complex eigenvalues, we say that the VT symmetry of H is broken. 

Lacking further information, one would expect (|2j) to have some complex roots. 
Thus, it was surprising to find that the class of PT-symmetric Hamiltonians ([1]) has an 
entirely real spectrum. In Refs. [2j [31 H] numerical evidence and detailed perturbative 
asymptotic analysis was presented to show that the eigenvalues of the Hamiltonian (OQ) 



Conjecture on the analyticity of VT -symmetric potentials 2 

are real when e > 0. A rigorous proof that these eigenvalues are all real was given by 
Dorey, Dunning, and Tateo 0, U\ ■ 

Using the WKB quantization condition 

" dxy/E - V(x) = (n + |) 7T, (3) 

'XI 

where x\ and X2 are the turning points [roots of V(x) — E = 0], one can understand 
heuristically why the eigenvalues of H in (TjQ) cease to be real when e < 0. The 
quantization condition (J3]) requires that there be a continuous integration contour from 
X\ to x 2 - Such a contour exists for e > 0, but as soon as e goes below 0, the contour 
joining x\ and x-i is broken by the existence of a branch cut in the complex- x plane and 
there is no longer a finite-length path connecting the turning points. 

The discovery that the eigenvalues of H in ([IT) are real led to a search for and 
subsequent study of other non-Hermitian PT-symmetric Hamiltonians whose spectra 
are also real [SI EH Efl El El EH E] • We emphasize that the reality of the eigenspectrum 
is an unusual property of a complex Hamiltonian and that many PT-symmetric 
Hamiltonians do not have entirely real spectra. For example, while the ix 2 y potential 
studied in Refs. [151 US] has a real ground-state energy, it has recently been found that 
some of the other eigenvalues are complex jTTJ . 

In this paper we conjecture that analyticity of the potential is a necessary (but 
not sufficient) criterion for a non-Dirac-Hermitian Hamiltonian to have real eigenvalues. 
This conjecture is based on extensive numerical studies in which we have found that, 
except in rare cases, a PT-symmetric Hamiltonian H = p 2 + V(x) does not have a real 
spectrum if its potential V(x) is not an analytic function of x. An example of such a 
nonanalytic PT-symmetric potential, which is discussed in Sec. 121 is 

V(x) = ix\x\. (4) 

We show in Sec. El that this potential has only one real eigenvalue. 

A heuristic explanation of the role played by analyticity can be based on the 
WKB contour integral in (j2j). For complex PT-symmetric potentials the derivation 
and application of this integral makes explicit use of the analyticity of the potential. 
Of course, the potential of a Hermitian Hamiltonian need not be analytic because its 
turning points lie on the real axis. Here, the integral for the WKB quantization condition 
is unambiguously taken along the real axis and does not need to be deformed into the 
complex plane. By contrast, the turning points for a complex potential are likely to be 
complex, and thus the contour for the WKB integral necessarily lies off the real axis. 
Giving up Hermiticity forces the quantization condition (J3J) into the complex plane and 
thus introduces the requirement of path independence and hence analyticity. 

The square- well potential studied by Znojil in Ref. [8] is a rare example of a complex 
nonanalytic "PT-symmetric potential having a real spectrum. This potential evades the 
conjecture above simply because there are no turning points at all; for the square-well 
potential there is no solution to the equation V(x) = E. 
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This paper is organized very simply: In Sec. [2] we examine the exactly solvable 
nonanalytic potential in (j3j) and in Sec. [3] we present numerical results for the class of 
nonanalytic PT-symmetric potentials 

V(x) = (ix) a \x\ b (a, b real). (5) 

In Sec. H] we make some brief concluding remarks. 

2. An exactly solvable nonanalytic potential 

In this section we consider the PT-symmetric Hamiltonian 

H = p 2 + ix\x\, (6) 

whose potential is a nonanalytic function of x. The Schrodinger eigenvalue differential 
equation associated with this Hamiltonian is 

d 2 \ 

-— + ix\x\-E\^{x) = 0, (7) 

where x is real and where the eigenfunction if)(x) is required to obey the boundary 
conditions that ip — > as i - >• ±oo. 

To solve this differential equation, we partition the real axis into two regions. In 
the region x > the differential equation (j7|) takes the form 

d 2 \ 

'— + ix 2 -E\^(x) = (8) 

and the exact solution is 

i>(t) = c x D v (xe i7T/8 Vf) + c 2 D u (-xe i7T/8 Vf) . (9) 

Here, D v is the parabolic cylinder function with 

v = \Ee-™ /A - |, (10) 

and C\ and c 2 are arbitrary constants. The boundary condition linx^+oo ip(x) = 
implies that c 2 = 0. Thus, for x > we have 

ip(x) = c x D u (xe i7T/8 V2) ■ (11) 

Similarly, in the region x < the differential equation (0) becomes 

d 2 \ 

■—-ix 2 -EU(x) = 0, (12) 

whose exact solution is 

V>(s) = dxDp (xe-^v^) + d 2 D ll {-xe" m,8 V2) . (13) 

Here, 

H = \e™ /A E - i (14) 
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and d\ and d 2 are arbitrary constants. The boundary condition lim, c _^_ 00 ■?/>(£) = 
implies that d\ = 0. Thus, for x < we have 

ij>{x) = d 2 D tx (-xe-**'*V2\ . (15) 

We must patch the two solutions (TTTT) and (TT5|) together at the origin x = 0. 
Continuity of ip(x) at x = implies that 

ciA,(0) = 4^(0), (16) 

and continuity of ip'(x) at x = implies that 

Cte^DUO) = -d 2 e-^ 8 D'^(0). (17) 

Taking the ratio of ( TT7I) and ( 1T6|) eliminates the constants ci and g?2 and gives the 
following exact equation for the eigenvalues: 

D v (0) D^O) ' [ ' 

This condition can be rewritten simply in terms of Gamma functions as 

e-/ 8 f T e . J + e— / 8 f r 6 . J = 0. (19) 

As required by PT symmetry, the secular equation ffT9l) is a real function of E. This is 
so because it is the sum of two terms that are complex conjugates of each other. 

To solve (TT9|) for E, we substitute E = Re E+i Im £? and take the real and imaginary 
parts of the resulting equation. We then plot in Fig. [T]the curves in the complex- E plane 
along which the real part of ffl9l) vanishes (solid line) and the imaginary part of ffl9l) 
vanishes (dotted line). [Of course, the condition of VT symmetry requires that the 
dotted line lie along the real-E axis. However, note that the real-E axis is not the only 
curve along which the imaginary part of (fl9|) vanishes.] 

The intersections of the solid and dotted lines in Fig. [1] are the eigenvalues of H 
in ([6]). Note that there is only one real eigenvalue; all other intersections occur in 
complex-conjugate pairs. The numerical value of the real eigenvalue is 

E = 1.258092 .... (20) 

3. Numerical study of a class of nonanalytic potentials 

In this section we examine numerically the eigenvalue differential equation for the 
complex PT-symmetric non-Hermitian Hamiltonian 

H = p 2 + [ix) a \x\\ (21) 

where a and h are real parameters. 

We begin by determining the appropriate PT-symmetric boundary conditions to 
be imposed on the eigenfunctions of H in (1211) . Using WKB analysis, we find the 
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Figure 1. Numerical solution to the secular equation (fill)) for the eigenvalues of the 
Hamiltonian H — p 2 +ix\x\. The solid line is the region in the complex- E plane where 
the real part of the secular determinant vanishes. The dotted line indicates where 
the imaginary part of the secular equation vanishes. The intersections of the solid 
and dotted lines are the eigenvalues. There is one real eigenvalue, which is located at 
E = 1.258092. All other eigenvalues come in complex-conjugate pairs. 



possible asymptotic behaviors of the eigenfunction ip(x) satisfying the time-independent 
Schrodinger equation 

<h ' ■■' V(x) =0. (22) 



dx 2 



+ (ix) 



\x\ b -E 



For example, when x > 0, the controlling factor of the asymptotic behavior of ip(x) as 

x — > +oo is 

2 



cxp 



_ i a/2 x (a+b+2)/2 



a + 6 + 2- - ■ < 23 > 

Thus, for a < 2 there exists a solution that grows exponentially and another that 
decays exponentially for large positive x. The same is true for large negative x so long 
as a < 2. To determine the eigenvalues for a < 2 we impose the boundary conditions 
that ip(x) — > as \x\ — > oo on the real-x axis. Note that because the potential is not 
an analytic function of x, the notion of Stokes' wedges in the complex-x plane in which 
the boundary conditions are imposed is not applicable. 

We have calculated the eigenvalues for various values of a and b and our results are 
listed in Tables [1] and [2] and plotted in Fig. [2j Clearly, when b is an even integer, the 
Hamiltonian in (1211) reduces to that in (JTJ) and it has an entirely real spectrum. However, 
for other values of b, when a/0 there are only a finite number of real eigenvalues. The 
number of real eigenvalues decreases as a increases, and increases as b increases. 
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Ei 


6 = 1/2 


6 = 1 


6 = 3/2 


6 = 2 


6 = 5/2 


6 = 3 


a = 


1.059617 


1.01879 


1.00118 


1 


1.00859 


1.02295 




1.833394 


2.33811 


2.70809 


3 


3.24223 


3.45056 




2.210015 


3.24820 


4.17714 


5 


5.72682 


6.37029 




2.550647 


4.08795 


5.58566 


7 


8.31328 


9.52208 




3.051182 


4.82010 


6.92282 


9 


10.9916 


12.8703 




3.253157 


5.52056 


8.22687 


11 


13.7342 


16.3694 




3.452132 


6.16331 


9.49059 


13 


16.5353 


20.0009 




3.623138 


7.37218 


10.7317 


15 


19.3837 


23.7455 




3.793400 


8.48849 


11.9453 


17 


22.2757 


27.5924 




3.943821 


9.53546 


13.1419 


19 


25.2052 


31.5308 


a = 1/2 


1.180777 


1.08693 


1.05583 


1.04896 


1.05404 


1.06568 




♦ 


3.19578 


3.27843 


3.43454 


3.59460 


3.74791 






4.4220 


5.36421 


6.05174 


6.64515 


7.17496 






♦ 


7.67568 


8.79101 


9.91884 


10.9735 








9.53919 


11.6207 


13.4256 


15.1112 








♦ 


14.5219 
17.4829 
20.4952 
23.5529 
26.6504 


17.0514 
20.8691 
24.7239 
28.8137 
32.7868 
37.2141 
41.0803 
46.2256 
♦ 


19.4889 
24.1139 
28.9111 
33.9218 
39.0482 
44.3936 
49.7770 
55.4476 
60.9963 
67.0561 
72.5848 
79.2958 
84.3126 
92.7345 
♦ 



Table 1. Real eigenvalues for the Hamiltonian H — p 2 + (ia;) a |a;| b in (f2"Tj) . where a = 
and a — 1/2. The symbol ♦ means that there are no more real eigenvalues, while . . . 
indicates that the spectrum is entirely real. 



4. Concluding remarks 



Most PT-symmetric potentials V(x) studied so far in the literature are special because 
they are analytic. In this paper we have explored a new class of nonanalytic VT- 



symmetric potentials of the form V(x) 



[ix ) \x\ 



and on the basis of numerical and 
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Ei 


6 = 1/2 


6 = 1 


6 = 3/2 


6 = 2 


6 = 5/2 


6 = 3 


a= 1 


1.446448 


1.25809 


1.18627 


1.15627 


1.14615 


1.14685 




♦ 


♦ 


4.21683 


4.10923 


4.13051 


4.19436 








6.93323 


7.56227 


7.95153 


8.30206 








♦ 


11.3144 
15.2916 
19.4515 
23.76667 
28.2175 
32.7891 
37.4698 


12.0844 
16.8072 
21.3065 
27.4779 
30.3268 
♦ 


12.9101 
18.1062 
23.5322 
29.6147 
35.3873 
42.9034 
47.4048 
♦ 


a = 3/2 


1.791941 


1.48873 


1.36338 


1.30151 


1.26993 


1.2550 




♦ 


♦ 


5.52801 


4.96979 


4.80096 


4.7494 








8.50818 


9.48003 


9.60759 


9.7042 








♦ 


14.5305 

19.9977 

25.8103 

31.9205 

38.2938 

44.904 

51.7304 


14.6672 
21.7069 

24.9567 


15.2406 

21.891 
28.1147 



Table 2. Same as in Table Q] except that a — 1 and a = 3/2. 

theoretical calculations we are led to conjecture that, except in rare cases, analyticity 
is an essential feature that is necessary for the Hamiltonian to have a real spectrum. 

We thank D. C. Brody for helpful discussions. CMB is supported by a grant from the U.S. Department 
of Energy. DWH receives financial support from Symplectic Ltd. 
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Figure 2. Real eigenvalues for the Hamiltonian H = p 2 + (ix) a \x\ b in (|2Tj) . The 
data for this figure is taken from Tables Q] and [2l The symbol ♦ indicates that there 
are no more real eigenvalues in the tower. Notice that the number of real eigenvalues 
increases with increasing b and decreases with increasing a. 
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